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Abstract 

Dispersion relations for fermions at high temperature and in a background magnetic 
field are calculated in two different ways. First from a straightforward one- loop cal- 
culation where, in the weak field limit, we find an expression closely related to the 
standard dispersion relations in the absence of the magnetic field. Secondly, we derive 
the dispersion relations directly from the Hard Thermal Loop effective action, which 
allows for an exact solution (i.e. to all orders in the external field), up to the last 
numerical integrals. 
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1 Introduction 



We know that thermal effects on dispersion relations are of extreme importance at high 
temperature, where the whole concept of propagating particles is drastically changed |T|, |2|. 
Since these effects are large and influence the mass shell conditions, it is also important 
to include them consistently to all orders and not only as first-order corrections. This is 
consistent to leading order in (eT) 2 . Both fermions and gauge bosons develop new branches 
and the zero-temperature masses become less important. New decay channels can open up 
because of the effective thermal masses of the propagating excitations. It is also interesting 
to study how strong external fields influence the physics of these effective modes. We have 
considered two ways of calculating the dispersion relations for fermions at high temperature 
in a constant background magnetic field. One method is to simply do the one-loop calcula- 
tion at finite temperature, using exact propagators in the external field Q . The other is to 
start from the already resummed effective action for Hard Thermal Loops 0, |5| (HTLs). It 
is an effective action that generates all HTLs at tree level and which contains the external 
fields to all orders. 



2 One-loop self-energy 

In Ref. [3] the fermionic dispersion relation at finite temperature in a constant magnetic 
field was calculated using the electron propagator in the Furry picture and Schwinger's 
proper-time method. In the high-temperature limit it turns out that Schwinger's proper- 
time formulation |J of the exact propagator is the easiest to use for the loop calculation. 
The thermal propagator can be constructed as j?| 

iSvac(p) - fF{po)[iS v3jC {p) - ^vacb) 



where for a magnetic field parallel to the z-direction[] 



r, / n f°° , e ieBsaz r / 2 taneSs 2 , 

»S»c(p) = J o ds^^ exp M [p, - - m + zs 




After performing the loop integrals in £(a/, x) = (x'|S|a;) it is possible to extract the gauge- 
invariant operator T,(p ,p z , U±), where = (p x — eA x ,p y — eA y ), in a gauge with zero 
^-component. (We shall use the gauge = (0,0, —Bx,0).) Then we take the weak field 
limit, but keep eB exactly wherever it is added linearly to the canonical momentum like 
U\-eBa z . The reason being that IT^ — eBa z = T/f±, so eB can be reabsorbed in which 
1 The sign convention here, that the particle has a positive charge e, differs from Ref. [3]. 
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need not be small. When the final expression for the self-energy is used in the effective 
Dirac equation, we obtain 



]/L-m-t(p ,p z ,U ± )}^ = 

s(p , n 2 )7op - r(p Q , H 2 )^ z p z - r(p , II 2 - eBa z )flj_ - m 
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The temperature dependence enters only through the thermal mass M. = e 2 T 2 /8. It is 
almost possible to guess the expression in Eq. (JJ) from the standard expression for the HTL 
Dirac equation Q. The usual momentum p M should be replaced with the gauge- invariant 
momentum ILj, but there is an ambiguity in replacing p 2 by IT 2 or by p. |/1 . The correct way 
follows from the calculations in Ref . [3] . 

In the massless limit (i.e. zero vacuum mass) the dispersion relation divides into a left- 
and a right-handed part. With a wave function \I> = (L, R) T in the chiral representation 
the two dispersion relations become 



L : {Es n +p z r n )(Es n ^i -p 2 r„_!) - 2eBnr 2 i 

2 

R : (Es n - p z r n )(Es n -! + p z r n -t) - 2eBnr 2 n _i 








(6) 
(7) 



where s n = s(E,p 2 , + eB(2n + 1)) and similarly for r n . These relations are valid for all n > 1, 
but in the lowest Landau level there is only one non-zero component for each of L and R, 
so the dispersion relations reduce to 



L : Es(E, p\ + eB) + p z r(E, p\ + eB) = 
R : Es{E,p 2 z + eB) - p z r(E,p 2 z + eB) =0 



(9) 



The numerical solution in the lowest Landau level is shown in Fig. [I] for the right-handed par- 
ticle. There are eight different branches corresponding to the states (L/R)x (particle/hole) 
x (positive/negative energy), for each given value of n > 1 and p y . For n = there are 
only half as many states, since only one spin orientation is possible in the lowest Landau 
level. Figure [l] only shows the right-handed branch, and the left-handed branch is obtained 
by reflection in p z . In fact, from Eqs. (|6|,|8|) we find the combined symmetries (L <-> R, 
Pz — Pz), (L <-> R,E <-> — E) and (E <-> —E,p z <-> — p z ). In particular, right-handed 
particles in the lowest Landau level can only propagate along the magnetic field (which 
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Figure 1: Dispersion relation and spectral weight for the right-handed branch in the lowest 
Landau level. All dimensionful parameters are given in units of the thermal mass M. 

points in the positive z-direction), while the left-handed particles and the right-handed 
holes propagate against the field. 

The spectral weight presented in Fig. [I] is defined by 



ui=Ei(p z ,n) 



(Tt [(X>(w,p*,7i)7o) 1 



-i 



(10) 



where D is the 4x4 matrix in Eq. (j3|) that remains after acting on an off-shell Landau 
level. For a given chirality, Zi(p z ) is obviously very asymmetric in p z since only particles 
propagate in one direction and the holes in the other. 



3 HTL effective action in a background field 

Instead of computing the one-loop self-energy as in Section 0, we can use the fact that the 
HTL effective action already contains the leading high-temperature contribution to all orders 
in the gauge field. The equation of motion that follows from a variation of the fermionic 
field should immediately give the dispersion relation. The HTL effective action for QED 
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can be written as 



C 



HTL 



4 



°-M*F pa 



u a u 



F t* 



(du) 2 

where the average (•) is defined by 



— u 



(f(u 0l u)} = J ^/(M) , 
where u is a spatial unit vector. The equation of motion for \l/ that follows is 



u-n 



* = o 



(ii) 



(12) 



(13) 



Equation (|H|) is a non-local and non-linear differential equation, which is, in general, very 
difficult to deal with. What makes this equation much less tractable than the thermal Dirac 
equation in the absence of the -B-field is that the average over u is difficult to perform 
because [IT M , Xl v ] = —ieF^ u ^ 0. Since the spatial symmetries of the system are unaltered 
by the thermal heat bath, we still expect the eigenfunctions to have the same spatial form as 
at zero temperature. In fact, after performing the w-integral in Eq. ( |i~3"|) the result can only 
be a function of the invariants 11^, p\ and p z , and the 7-structure has to be proportional to 
7II_l, 70P0 an d r y z p z - We shall therefore compute the matrix elements 



(3V 



u ■ n 



l*«> 



between the vacuum eigenstates 
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:py 



(14) 
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(16) 



where k = {po,n,p y ,p z } and i/" n [x] are Hermite polynomials. These states form a complete 
set of functions in four dimensions when the energy is off shell. In the chiral representa- 
tion suitable spinors can be formed from $ K as \& K = diag[$ K , $ K , where \ is 
an undetermined space-time-independent spinor. The vacuum Dirac operator in Eq. ( |13|) 
obviously gives an eigenvalue when acting on but it is more difficult to determine the 
action of the thermal part since <3> K cannot be an eigenfunction to u • II for all u. One way to 
calculate the matrix element in Eq. (|T4|) is to find a basis such that v ■ TL\v p ) = v ■ p\v p ) and 
insert a unit operator / d 4 p\v p ) (v p \. This unit operator is, of course, independent of v after 
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the p-integration, so in particular we can choose v = u and change the order of integrations. 
After computing the matrix elements in Eq. fll4]) we find indeed that they are diagonal in k 
for uo and u z , and have a mixing with the first subdiagonals for u x and u y . Define (uo )Z) ±) 
by 

(^'l(^l)l^) = (27r) 3 <WK*>« , (17) 
= (2n) 3 5 K ^ T i(u ± ) K , (18) 

and k =F 1 = {po, l,p y ,p z }. These are exactly the components that occur naturally when 
we include the 7-matrices in the chiral representation. The calculation of (uo, z ,±) is a bit 
lengthy but straightforward, and the result reads 



n!V2vr 
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(u z ) K 



\= f°° dsH 2 n (s)e~ s2 / 2 
. 2ix J —00 



n\y27c J-00 

f Pz Po(2p 2 z ~ eBs 2 ) Po+p z 



I P 

p z (2E 2 



p 2 ) s 



V2eB 



2p 4 



Po~P 2 



arctan • 



Po-Pz 



2Jp 2 o~P 2 



(20) 



i f°° 

(u+) K = . = / dsH n (s)H n - 1 (s)e' 

x 2im\(n - 1)\ J ~°° 



■s 2 /2 



\s\/2eB Esp z V2eB , p + p z 

x \ o a l n 

I 2p l 2p A p - p z 



+ 2 ^M-f -^'arctan "f* ) , (21) 
2p 4 JpI~p 2 2yPo-p 2 



(«-)« = -<u+)«+i , (22) 

where p 2 = p 2 . + eBs 2 /2. With these definitions the Dirac equation effectively reduces to a 
4x4 matrix in the spinor indices, since the other quantum numbers have been diagonalized. 
In the massless limit we take the determinant to find the dispersion relations. They factorize 
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again as in Eqs. and the equations for the right-handed component are 

n>l: (po -p z - M 2 {(u ) K - (u z ) K f) 

x (po + p z - M 2 {(u ) K -i + (u z ) K -i)J 

\ 2 ^ ' 

- [V2eBn -iM 2 (u + ) K ) = , 
n = 0: Po-Pz- M 2 {(u ) K - (u z ) K ) = . 

In general there are imaginary parts in the functions (uq >Z) ±) k , which make the spectral 
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Figure 2: Comparison of the dispersion relation from the HTL effective action and the weak 
field approximation. All dimensionful parameters are given in units of the thermal mass Ai. 

functions more complicated. The imaginary parts of Eqs. ([H]) to fl22|) are determined by 
the analytic continuation p — > p + ie for positive po; this amounts to the replacement 



1 



pl-p 2 



arctan 



sV2c~B 



In 



sV^eB + 2 Jp 2 - p 2 



VK 



2Jp 2 -p 2 sV2eB-2Jp*- P 2 



P 2 ~Po 



(24) 



6 



for p 2 > Pq. It is anyway useful to solve Eq. (|2"3|), ignoring the imaginary part, since the real 
part indicates where the spectral functions are peaked, at least if the imaginary parts are 
small enough. This can conveniently be done numerically as all the integrals in Eqs. ( |19|) 
to ( |2"2"D are well convergent. The result for the lowest Landau level in a weak magnetic field 
(eB = 0.2 A4 2 ) is shown in Fig. |2|, and it agrees rather well with the weak-field result from 
Eq. @. 
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Figure 3: Dispersion relation for the right-handed branch in the lowest Landau level (n = 0). 
As the B-Eeld increases thermal effects become less important and the dispersion approaches the 
light cone, which is indicated by solid lines. All dimensionful parameters are given in units of the 
thermal mass Ai . 

It is, of course, more interesting to see what happens at larger field strength, which 
cannot be treated by Eq. (U). The dispersion relations for several field strengths are shown 
in Fig. Apart from the changes in the particle branch there is a new branch coming 
from the light cone, which eventually joins the hole branch and disappears. This new 
branch can be understood mathematically, from Eq. (|23|) , by studying it close to the light 
cone. In the absence of the 5-field, the hole branch exists because the logarithm in Eq. (f|) 
becomes dominant close to the light cone, and the sign is such that it allows for a positive 
energy solution to the part of that Dirac equation which normally only gives the antiparticle 
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solution. In the present case, there is a compensation from the new terms in Eq. (^3]), which 
change the behaviour close to the light cone and a new branch can exist. As the i?-field 
increases the hole branch and its partner become less extended, and the particle branch 
approaches a massless mode. This is physically very reasonable since for very strong field 
strengths the thermal effects should disappear. 

At this point it is worth discussing the approximations involved. The HTL effective 
action is derived under the assumption that the temperature is much larger than the external 
momenta. Here, the momentum is at least II ~ \fe~B, which should be kept smaller than T. 
But, if we consider the coupling constant e to be very small this approximation should be 
valid even for eB ~ M 2 . On the other hand, when II 3> T the vacuum part is dominant and 
the dispersion relation is still approximately valid. The other approximation was to neglect 
the imaginary part. This approximation becomes worse for increasing eB and cannot be 
motivated for large eB. A new massless excitation, such as the new hole partner, could have 
important physical consequences only if its spectral weight is non-negligible. This remains 
to be studied and it can only be done correctly using the full spectral function. 
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